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Abstract: 

Spectroscopy of planetary nebulae (PNe) provides the means to investigate s-process enrichments 
of neutron(n)-capture elements that cannot be detected in asymptotic giant branch (AGB) stars. 
However, accurate abundance determinations of these elements present a challenge. Corrections for 
unobserved ions can be large and uncertain, since in many PNe only one ion of a given n-capture 
element has been detected. Furthermore, the atomic data governing the ionization balance of these 
species are not well-determined, inhibiting the derivation of accurate ionization corrections. We present 
initial results of a program that addresses these challenges. Deep high resolution optical spectroscopy 
of ~20 PNe has been performed to detect emission lines from trans-iron species including Se, Br, Kr, 
Rb, and Xe. The optical spectral region provides access to multiple ions of these elements, which 
reduces the magnitude and importance of uncertainties in the ionization corrections. In addition, 
experimental and theoretical efforts are providing determinations of the photoionization cross-sections 
and recombination rate coefficients of Se, Kr, and Xe ions. These new atomic data will make it possible 
to derive robust ionization corrections for these elements. Together, our observational and atomic data 
results will enable n-capture element abundances to be determined with unprecedented accuracy in 
ionized nebulae. 
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1 Introduction 

Approximately half of the neutron(n)-capture elements 
(atomic number Z > 30) in the Universe are created 
by slow n-capture nucleosynthesis (the "s-process"). 
The s-process can occur in low- and intermediate- mass 
stars (1-8 Mq), the progenitors of planetary nebulae 
(PNe), during the thermally-pulsing asymptotic giant 
branch (AGB) phase. Free neutrons are released by 
the reaction 13 C(a,n) ls O — or 22 Ne(a,n) 25 Mg in more 
massive AGB stars (> 3.5 Mq) — in the intershell re- 
gion between the H- and He-burning shells. Fe-peak 
nuclei experience a series of n-captures interspersed 
with /3-decays to transform into isotopes of heavier el- 
ements. The enriched material is transported to the 



stellar envelope via convective dredge-up, and is ex- 
pelled into the ambient interstellar medium by stellar 
winds and PN ejection ijBusso et al. IH999T ). 

Nebular spectroscopy uniquely reveals information 
about n-capture nucleosynthesis that cannot be ob- 
tained from stellar spectra. For example, the light- 
est n-capture elements (Z = 30-36) and noble gases 
are not well-studied in their sites of synthesis, due 
to the difficulty in detecting these species in evolved 
stars or supernova remnants (e.g., Wall erstein et al~l 
1995). In contrast, these elements are readily detected 
in PNe. Spectroscopy of PNe also enables investiga- 
tions of s-process nucleosynthesis in intermediate- mass 
stars (M = 3.5-8 Mq). These stars are obscured 
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by dusty, optically thick circumstellar envelopes dur- 
ing the AGB and post-AGB phases (jHabing I I"l996l : 
IVan Winckel Il2003h that preclude optical or UV spec- 
troscopy. Hence, little is known about their contribu- 
tion t o s-process nuclei in the Universe l|Travaglio et al. 1 
1 20041 ) . These stars produce Type I PNe (cha racterized 
by He and N enrichments; iPeimbertl Il978f ) in which 
n-capture elements can be detected. 

The low cosmic a bundances of n-capture elements 
(jAsplund et al. Il2005|) prevented their identification in 
nebular spectra until lPeauignot fc Baluteau I (|1994l . here- 
after PB94) discovered emission lines of Se, Br, Kr, Xe, 
and possibly other trans-iron elements in the optical 
spectrum of the PN NGC 7027. Many of these iden- 
tifications were later confir med with highe r resol ution 
spectroscopy perform ed by IZhang et al. I <|2005T) and 
ISharpee et al. I (120071 ). PB94's analysis led lDinerstein I 
(|200ll ) to identify two anonymous emission lines in the 
K band spectra of PNe as fine-structur e transitions of 
[Se IV ] and [Kr III]. Shortly thereafter. Isterling et al~l 
(2002) identified Ge III A1088.46 in absorption in the 
UV spectra of four PNe, and found that these objects 
exhibit Ge enrichments indicative of s-process nucle- 
osynthesis in their progenitor stars. Together, these 
studies determined abundances of n-capture elements 
in only 11 PNe, too small a sample to characterize 
s-process enrichments in PNe as a population. 

ISterling fc Dinerstein I [|200Sl . hereafter SD08) per- 
formed the first large-scale survey of n-capture ele- 
ments in PNe by detecting the near-infrared emission 
lines [Kr III] 2.199 and/or [Se IV] 2.287 /im in 81 
of 120 Galactic PNe. To correct for the abundances 
of unobserved Se and Kr ions, analytical ionization 
correction factors (ICFs) were derived from a larg e 
grid of photoionization models (jSterling et al. | |2007). 
We derived Se and Kr abundances (or upper limits) 
in each object, increasing the number of PNe with 
known n-capture element abundances by nearly an or- 
der of magnitude. Se and Kr were found to be en- 
riched by a factor of two or more in ~45% of the ob- 
served PNe, indicating efficient s-process nucleosyn- 
thesis and convective dredge-up in their progenitor 
stars. Both exhibit a positive correlation with the 
C/O ratio, as ha s been found for other n -capture el- 
ements in AGB (ISmith fc Lambert1ll990l ) and post- 
AGB stars (| Van Winckel 1 120031 ) . Type I PNe, which 
have intermediate-mass progenitors, display little or no 
enrichment of Se and Kr compared to other PN sub- 
classes. These trends have be en shown to be consis tent 
with theoretical predictions (jKarakas et al. II2009T ). 

The K band Se and Kr lines were recently detected 
i n thr ee nearby extragalactic PNe bv lDinerstein et al. I 
( 2009 ) , providing insight to the nucleosynthesis of these 
elements in low-metallicity environments for the first 
time. The PN Hen 2-436 in the Sagittarius dwarf 
spheroidal galaxy exhibits large enrichments of Se and 
Kr (Figure 1), comparable to the most highly enriched 
Galactic PNe observed by SD08, while two PNe ob- 
served in the Magellanic Clouds do not appear to be 
enriched. 

The Se and Kr abundances derived in these stud- 
ies are typically accurate to within a f actor of two or 
three, based on the error analyses of Ste rling et al. I 
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Figure 1: K band spectrum of the PN Hen 2-436 in 
the Sagittarius dwarf spheroidal galaxy, measured 
with the GNIRS spectrogra ph on Gemini South 
(from lDinerstein et al. |[2009D . Both [Kr III] 2.199 
and [Se IV] 2.287 ^m were detected, and are iden- 
tified in the figure. 



(120071 ) and SD08, and similar uncertainties can be ex- 
pected for other n-capture elements observed in ion- 
ized nebulae. These uncertainties stem from (1) the 
detection of only one ion each of Se and Kr, which 
can lead to large and uncertain ICFs; and (2) the 
poorly known atomic data that controls the ioniza- 
tion balance of these elements, which precludes the 
derivation of accurate ICFs. We present initial results 
from a multi-disciplinary investigation, comprised of 
new optical observations and atomic data determina- 
tions, designed to address these sources of uncertainty 
and enable much more accurate abundance analyses of 
n-capture elements in ionized nebulae. 

2 Optical Observations 

Optical spectroscopy provides access to transitions from 
several ions of Br, Kr, Rb, Xe, and (when combined 
with near-infrared data) Se. Detecting multiple ions of 
these elements reduces the magnitude and importance 
of uncertainties in their ICFs, and hence in the derived 
elemental abundances. Moreover, Xe is a "heavy-s" 
element that lies near the s-process peak at Ba. The 
abundance of Xe relative to those of lighter n-capture 
elements can be used to determine the neutron ex- 
posure (i.e., the time-averaged neutron flux) experi- 
enced by Fe-peak nuclei during s-process nucleosynthe- 
sis. At a given metallicity, the neutron exposure con- 
trols the ele ment-by-element distribut ion of s-process 
enrichments j|Busso et al. II 19991 . [200lT ). 

The robust identification of optical n-capture el- 
ement lines in nebular spectra requires high spectral 
resolution in order to distinguish them from nearby 
transitions of more abundant elements and/or telluric 
features. We are conducting optical observations of s- 
process enriched Galactic PNe selected from the near- 
infrared sample of SD08, using the 2D Coude cross- 
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Figure 2: Selected Kr and Xe emission lines detected in the PN J 900 (top panels), and Br and Rb lines 
in NGC 6644 (lower left panels) and NGC 3242 (lower right panels). Telluric emission lines have not 
been removed, and appear as narrow features within some of the panels. Note that wavelengths are not 
corrected for systemic nebular velocities. 



dispersed echelle spectrometer (jTull et al. I Il995l ) on 
the 2.7-m Harlan J. Smith Telescope at McDonald 
Observatory. Using two grating tilt settings and a 
1.8" -wide slit, this instrument provides nearly com- 
plete spectral coverage in the wavelength range 3 700- 
10 400 A, at a resolution R — 36 700 (sufficiently large 
to resolve many n-capture element lines from potential 
blends). 

To date, we have completed optical observations 
of approximately 17 PNe. These objects display rich 
spectra, with ~100 emission lines in the sparsest spec- 
trum to over 400 lines in NGC 3242. Kr emission has 
been detected in all but one of these objects, and sev- 
eral exhibit lines from multiple Kr ions. At least one 
Xe emission line has been detected in about half of 
these PNe, and Se, Br, and Rb lines have been tenta- 
tively identified in 3-4 objects each. Figure 2 displays 
selected Kr, Xe, Br, and Rb emission lines in targets 
from our sample. 

3 Atomic Data Calculations 

Abundance determinations of n-capture elements in 
ionized nebulae are affected by uncertainties in the 
atomic data governing the ionization balance of these 
species — in particular, photoionization (PI) cross- 
sections and rate coefficients for radiative recombina- 
tion (RR) , dielectronic recombination (DR) , and charge 



transfer (CT). These data are needed to derive ac- 
curate ICFs to correct for the abundances of unob- 
served ions. Atomic data uncertainties can lead to 
Se and Kr abundanc e uncertainties of a factor of two 
([Sterling et al. Il2007l ) . and similar uncertainties can be 
expected for other n-capture elements. 

Wit h the atomic s tructure code AUTOSTRUC- 
TURE (jBadnell iri986t ). we are computing PI cross- 
sections and RR and DR rate coefficients for the first 
five ions of Se, Kr, and Xe, the t hree most widely de- 
tected n-capture elements i n PNe (Pequienot & Bal uteau I 
119941 : ISharpee et al. I [20071 . SD08). This code can be 
used to efficiently compute multi-configuration distorted 
wave PI cross-sections and radiative and autoioniza- 
tion rates. Relativistic effects are accounted for via 
the Breit-Pauli formalism and semi-relativistic radial 
functions. 

We are computing the electronic structure of each 
ion using multiple configurations within the intermedi- 
ate coupling scheme. These configuration expansions 
are utilized to compute PI cross-sections out of the 
valence shells, which are resolved in both initial and 
final state and can be combined to determine the to- 
tal cross-section of each ion. In Figure 3, we show 
the computed Xe 3+ PI cross-section near the ground 
state threshold regi on, plotted against the experimen- 
tal measurements o f lBizau et al"Tl|2006l ). RR is related 
to the non-resonant portion of the PI cross-sections, 
and can be computed from the principal of detailed 
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Figure 3: The calculated Xe 3+ PI cross-section 
(solid line), using a 5 configuration expansion for 
Xe 3+ and 6 configurations for Xe 4+ . Note that 
the resonance contribution has not been added. 
For comparison, we show th e experimental Xe 3+ 
PI cross-section measured bv lBizau et al. I ((2006). 



balance. Similarly, DR can be related to the reso- 
nant portion. Our DR rate coefficients are affected by 
uncert ainties in the e nergies of low-lying autoionizing 
levels ijFerland I l2003h . which have not been spectro- 
scopically determined for these ions. However, these 
uncertainties can be alleviated with our experimental 
PI cross-section measurements (§4). 

For each of these processes, we are estimating the 
uncertainties in our computations in order to evaluate 
their effect on elemental abundance derivations (§5). 
For example, uncertainties in the PI cross-sections and 
RR rate coefficients can be estimated by altering the 
configuration expansions utilized for the initial and 
target ions. DR uncertainties are estimated by shifting 
the continuum level relative to the near-threshold en- 
ergy levels, thereby varying the energies and number 
of autoionizing channels available for DR. 



Experimental Photoionization 
Cross-Section Measurements 
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Figure 4: Top panel: The Se + photoionization 
cross -section measured at the ALS ([Sterling et al. I 
2009). The circles with error bars represent abso- 
lute measurements performed at discrete photon 
energies. Bottom panel: Magnified view of the 
Se + PI cross-section in the metastable region (the 
ground state ionization threshold is at 21.2 eV), 
with resonances and Rydberg series identified. 



To benchmark our atomic data calculations, we are 
performing experimental absolute PI cross-section mea- 
surements of Se, Kr, and Xe ions. These data con- 
strain both our PI and DR calculations, since PI res- 
onances near the ionization threshold correspond to 
low-temperature DR resonances. Measuring the en- 
ergies of these resonances allows us to normalize the 
unknown energies of autoionization resonances in our 
DR calculations. 

These measurements are being conducted at the 
Advanced Light Source (ALS) synchrotron radiation 
facility at the Lawrence Berkeley National Labora- 
tory, with the ion- photon merged beams endstation on 
Beamline 10.0.1.2 ([Covington et al. l2002l ; lAguilar et al. I 
2003). In these measurements, a beam of atomic ions 



is merged with a high-energy photon beam to produce 
photoions. EUV and soft X-ray photons are produced 
by electrons accelerated to 1.9 GeV that pass through 
an undulator located in a straight section of the storage 
ring lattice. At the beamline, the photon beam energy 
and resolution are controlled by a spherical-grating 
monochromator with three interchangeable gratings and 
adjustable slits. Atomic ions are produced in an electron- 
cyclotron resonance source, and accelerated by a po- 
tential (typically 6 keV). The desired charge state is 
selected with an analyzing magnet, and the ion beam 
is merged with the counter-propagating photon beam. 
Photoions are produced along the merged beam path, 
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Figure 5: The Se 3+ photoionization cross-section 
measured at the ALS. Circles with error bars rep- 
resent absolute measurements performed at dis- 
crete photon energies. 



and are directed to a detector by a demerging magnet. 

This apparatus has two operational modes: spec- 
troscopy and absolute mode. For spectroscopy, pho- 
toion yields produced along the entire merged path 
are recorded in arbitrary units as a function of pho- 
ton energy. To place these yields in physical units 
(absolute mode), a potential is placed on a cylindri- 
cal "interaction region," which has a precisely defined 
length. This potential energy-tags all photoions pro- 
duced within the interaction region, and the demerg- 
ing magnet setting is altered so that only the energy- 
tagged ions are detected. By determining the overlap 
of the photon and ion beams in the interaction region, 
absolute PI cross-sections can be measured at discrete 
energies. These measurements are used to normalize 
the photoion yields onto an absolute scale. 

We have measured absolute PI cross-sections near 
the ionization thresholds of Se + ( Ster ling et al. I 2009, 
Figure 4), Se 2+ , Se 3+ (Figure 5), Xe 2+ , and the pho- 
toion yield of Xe + . PI cross-sections of K r ions hav e al- 
ready been measured l|Lu et al. ll2006al lbl:lLu 2006), as 
have those of more highly charged Xe ions (|Bizau et al~l 
2006). Unfortunately, the Se 4+ PI cross-section can- 
not be measured with this apparatus due to an Auger 
resonance effect (also seen in the isoelectronic Kr 6+ ). 

The analysis of these measurements is complicated 
by the presence of meta-stable states in the primary 
ion beam, which may contribute to the measured cross- 
sections along with ground state ions. It is neces- 
sary to measure the PI cross-sections at energies below 
the ground state threshold (i.e., down to those of the 
metastable states) in order to estimate the fractional 
content of these states in the ion beam. 

The measured cross-sections are interpreted with 
the aid of the fully relativistic 7?-matrix code DARCQ 



(jNorrington fc Grant Ill987l : IWiiesundera et al. Ill99lh 
in order to identify resonance features and to disen- 
tangle the contributions of the ground and metastable 
states (e.g., Figure 6). A parallel ve rsion of the PARC 
suite o f codes has been developed bv lBallance fc G riffin 
(2006), and was recently modified to efficiently com- 
pute numerous bound- free dipole matrix elements, dis- 
tributed over an equivalent number of processors to 
enable photoionization calcul ations to be carried out 
at a high degree of accuracy (|Ballance et al. II2008T ). 



5 Implementation 

The new atomic data will be incorporated into the 
atomic database of t he photoionization code Cloudy 
dFerland et al. 



1998). Following the methodology of 



20071 ) . we will compute a grid of mod- 



ISterling et al. 
els to determine analytical ICFs for Se, Kr, and Xe. 
In addition, we will use Cloudy to test the effects that 
our estimated atomic data uncertainties have on the 
derived elemental abundances. This quantitative anal- 
ysis will clarify the ionic systems and atomic processes 
that require further attention in future investigations. 
The ICFs will be applied to our optical (§2) and near- 
infrared spectra (SD08) to derive accurate n-capture 
element abundances in PNe, in order to study the de- 
tails of s-process nucleosynthesis and convective mix- 
ing in PN progenitor stars. Combined, our optical ob- 
servations and atomic data determinations will allow 
us to determine n-capture element abundances and s- 
process enrichment factors in ionized nebulae with un- 
precedented accuracy. 
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